Study of the P-wave charmonium state XcJ in ^(25') decays 
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Abstract 



The processes ip(2S) — > 77r + 7r~, ^K + K~ and 7pp have been studied using 
a sample of 3.79 x 10 6 ip(2S) decays. We determine the total width of the 
XcO to be r^°' = 14.3 ± 2.0 ± 3.0 MeV. We present the first measurement of 
the branching fraction B(x c o —> Pp) = (15.9 ± 4.3 ± 5.3) x 10~ 5 , where the 
first error is statistical and the second one systematic. Branching fractions of 
Xc0,2 — ► vr + 7r _ and K + K~ are also reported. 
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Typeset using REVTgX 



The hadronic decay rates of the P-wave quarkonium states provide tests of perturbative 
QCD. Recently, a systematic approach to the treatment of the infrared ambiguities in the 
calculation of the production and decays of these states has been developed [|TJ. However, 
existing experimental information on the triplet P-wave cc states (xco,i>2)> especially the 
J = Xco, is not adequate for testing the predictions of this new theory 

In particular, the total width of the Xco is a quantity of considerable interest. The two 
existing measurements have large errors and only marginal consistency ||. Also of interest is 
the decay Xco PP, which is forbidden in the limit of massless helicity conservation |J and 
has been calculated by many different models ||,|5| . Here the only existing measurement is an 
upper limit on the partial width that does not seriously constrain the theory || . Calculations 
of the branching fractions for other exclusive Xcj decays, such as Xcj — ► 0, has 

revealed orders-of-magnitude discrepancies with the data reported by early experiments. For 
these reasons, measurements of these properties of the Xcj states with improved precision 
are very useful. 

In this paper we report a measurement of determined from an analysis of exclusive 
ip(2S) — > 77r + 7r~ and , -fK + K~ decays seen in the Beijing Spectrometer (BES) at the Beijing 
Electron Positron Collider (BEPC). We also report a first measurement of the Xco —> PP 
branching fraction, and improved precision on the branching fractions for Xco,2 vr + 7r~ and 
K+K~. 

The BES is a conventional solenoidal magnet detector that is described in detail in 
Ref. ||. A four-layer central drift chamber (CDC) surrounding the beampipe provides 
trigger information. A forty-layer main drift chamber (MDC), located radially outside the 
CDC, provides trajectory and energy loss (dE/dx) information for charged tracks over 85% 
of the total solid angle. The momentum resolution is <r p /p = 0.017a/1 + p' 2 (p in GeV/c), 
and the dE/dx resolution for hadron tracks is ~ 11%. An array of 48 scintillation counters 
surrounding the MDC measures the time of flight (TOF) of charged tracks with a resolution 
of ~ 450 ps for hadrons. Radially outside the TOF system is a 12 radiation length, lead- 
gas barrel shower counter (BSC). This measures the energies of electrons and photons over 
~ 80% of the total solid angle with an energy resolution of ue/E = 22%/\fE (E in GeV). 
Outside of the solenoidal coil, which provides a 0.4 T magnetic field over the tracking volume, 
is an iron flux return that is instrumented with three double layers of counters that identify 
muons of momentum greater than 0.5 GeV/c. 

We study Xc states produced by the reaction e + e~ — > i/j(2S) — > jXc in a data sample 
corresponding to a total of 3.79 x 10 6 ip{2S) decays ||. For the r^ determination reported 
here we use the paired pseudoscalar meson decay modes of Xco,2 —> vr + 7r~ and K + K~ . Using 
the particle identification capabilities of the detector and four-constraint kinematic fits, we 
can get relatively pure event samples. Moreover, since the decays of the Xci and the rf c 
to 7r + 7r~ or K + K~ are forbidden by parity conservation, the Xco and Xc2 signals in these 
channels are free of distortions due to possible contamination of these other states. The 
effects of cross-contamination between the it + tt~ and K + K~ event samples are estimated 
by Monte Carlo (MC) simulations and corrected accordingly. The total width of the Xc2 has 
been precisely measured to be r^°* = 2.00 ± 0.18 MeV [plj, which is much narrower than 
our experimental resolution at M Xc2 of 7.83 MeV. Thus, the strong Xc2 — > k + -k~ signal in 
our data is used to provide a direct experimental determination of our resolution. We only 
rely on the MC simulation to determine how the resolution changes between M Xc2 and M Xc0 . 
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We select ip{2S) — > 77r + 7r~, , -fK + K~ and 7pp by imposing the following selection criteria. 

A cluster of deposited energy in the BSC is regarded as a photon candidate if: (1) the 
angle between the nearest charged track and the cluster in the r<p plane is greater than 
15°; (2) the energy of the cluster is greater than 20 MeV and some energy is deposited in 
the first 6 radiation lengths of the counter; and (3) the angle determined from the cluster 
development in the BSC agrees with that determined from the relative position of the shower 
location and the interaction point to within 37°. At least one and at most three photon 
candidates are allowed in an event. The candidate with the largest BSC energy is assumed 
to be the photon radiated from the ip{2S). 

In addition, we require that the event has two oppositely signed charged tracks in the 
MDC that both have at least 13 good hits and are well fit to a three-dimensional helix. 
Events with tracks where the dE/dx measured in the MDC and the shower properties in 
the BSC are consistent with electrons are rejected. For each track, the TOF and dE/dx 
measurements are used to assign probabilities that the particle is a pion, kaon and proton 
(Prob n , Probx, Prob p ). We require both tracks to have Prob n > 0.01 (for 7r + 7r~) or Probx > 
0.01 (for K + K~) or Prob p > 0.05 (for pp). In addition we do four-constraint kinematic fits 
to the hypotheses ip(2S) — ► •yn + 7r~ , ip(2S) — > r yK + K~ and ip(2S) — > jpp and require the x 2 
probability of the fit P x 2 to be greater than 0.01 for 7r + 7r~ or K + K~ and greater than 0.05 
for pp. 

There is some background from J/ip — > /i + yU~, where the J/ip is produced by cascade 
ip{2S) to J/ip decays. To reduce this, we reject events where the response of the muon 
detection system is consistent with the two charged tracks being muons. The surviving 
fi + fi~ background events do not populate the 7r + 7r~ invariant mass distribution near the Xco 
or Xc2 masses. In the K + K~ mass distribution, however, they populate the region in the 
lower mass side of Xco, an d cause an abnormal distribution. For the pp sample, the 
background level is significant. For this channel, in order to insure that both tracks are 
directed at the sensitive area of the muon detection system, we require | cos^tjqI < 0.65 
for both the p and the p track. 

To distinguish 77r + 7r _ from , -fK + K~ , we define 

Pr obau P ~ = Prob(xl\\, nd 4ll)' 

where x\\\ = xlc + Xtof + X 2 dE /dx and nd f&\\ = nd hc + nd froF + nd f dE/dx are the total 
X 2 and the corresponding number of degrees of freedom of the x 2 distribution. Here x\c-> 
Xtof an d XdE/dx correspond to the x 2 values from the 4-constraint kinematic fit, the TOF 
measurements for the 7r or K hypothesis, and the dE/dx measurements for the 7r or K 
hypothesis, respectively, and ndf 4C , ndf TOF and ndf dE j dx are the corresponding numbers of 
degrees of freedom. If Prob a ^ l 7T > Prob^ K , the event is categorized as a ^ti + tt~ event, 
and if Probst K > Pfob^^ , it is categorized as a r )K + K~ event. 

Figures [3] and |] show the 7r + 7r~ and K + K~ invariant mass distributions after the im- 
position of all the above-listed selection requirements. In these plots, the mass values cor- 
responding to the Xco an d Xc2 peaks are lower in the 7r + 7r~ channel and higher for K + K~, 
indicating the presence of some remaining cross contamination between the two samples, 
which must be accounted for in the determination of the Xco parameters. From Fig. ^, it 
is apparent that the Xc2 —> K + K~ sample is statistically limited. We therefore use only 
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the Xc2 —* 7r + 7r _ signal to calibrate the mass resolution. In Fig. |3], the pp invariant mass 
distribution, there is a clear Xco signal and evidence for the Xci and Xc2- 

We use Monte Carlo simulated data to determine the n + ir~ and K + K~ cross contami- 
nation probabilities, the detection efficiencies and the mass resolutions. We generate events 
assuming that the reaction ip(2S) — > JXcJ is a pure El transition. The decays Xcj — > pseu- 
doscalar meson pairs and Xco —> PP have only one independent helicity amplitude and are 
thus unambiguous O]. For x&,2 ~~ > VP decays, there are no available experimental data on 
the helicity amplitudes, and we use an isotropic distribution. (Our pp event samples are too 
small to permit a helicity amplitude analysis.) 

We subject the MC-generated events to the same selection process as is used for the 
data and determine the detection efficiencies for each mode. For the pp mode, the detection 
efficiencies and the error caused by the limited statistics of the Monte Carlo sample are 
e Xc0 = (27.1 ± 0.6)%, e Xcl = (30.3 ± 0.7)% and e Xc2 = (27.6 ± 0.6)%, and mass resolutions 
at the Xco, Xci and Xci equal to 7.3, 6.8 and 6.7 MeV, respectively. For the 7r + 7r~ and 
K + K~ modes, the simulation shows that the mass dependence of the detection efficiency 
is small and the mass resolution function is very nearly Gaussian. We compensate for 
the distortion of the mass spectra due to the the tt + tt~-K + K~ cross contamination by 
calibrating the mass resolution derived from the MC simulation with the Xc2 ~~ > 7r + 7r~ 
line shape seen with the data. The efficiencies are e Xc0 ^ n + n - = (36.9 ± 0.3)%, £ Xc2 _ +7r + 7r - = 
(38.9±0.5)%, e Xc0 _> K+K - = (32.8±0.3)% and e Xc2 ^ K + K - = (34.9±0.5)%, and the probability 
for Xco {Xci) — ¥ K + K~ events to be categorized as n + n^ is (5.6 ± 0.2)% ((6.0 ± 0.2)%), and 
that for Xco {Xci) -> events to be selected as K + K~ is (7.1 ± 0.2)% ((7.4 ± 0.3)%), 

where the error is from the statistics of the Monte Carlo sample. 

The invariant mass distributions in Figs. [I], [| and |^ are fit by using an unbinned maxi- 
mum likelihood algorithm. For the pp channel, the invariant mass region between 3.20 and 
3.64 GeV is fit with three Breit-Wigner resonances plus a linear background function. The 
Breit-Wigner resonance width for the Xco is fixed at 14.3 MeV, the value determined from 
an analysis of Xco — ► 7r + 7r~ decays described below; those for the Xci and Xc2 are fixed at the 



PDG values ||10|| . The resonances are smeared by Gaussian functions with rms widths fixed 
at the MC-determined mass resolution values. The fit result, shown as the curve in Fig. |3|, 
gives 15.2 ± 4.1, 4.2 ± 2.2 and 4.7 ± 2.5 events for the Xco, Xci and Xc2 states, respectively. 

For the 7r + 7r~ channel, we first fit the invariant mass region between 3.5 and 3.6 GeV 
with a Breit-Wigner resonance with T Xc2 fixed at the PDG value of 2.00 MeV, smeared by 
a Gaussian resolution function with an rms width that is allowed to float. We also include 
a linear background function. The fit results in a Xc2 mass resolution of 7.83 ± 1.04 MeV, 
which is slightly higher than the MC result of 6.31 ± 0.11 MeV. We scale the MC value for 
the mass resolution at the Xco (8.12 ± 0.23 MeV) by the ratio of the fitted MC results at the 
Xc2 and get a mass resolution at the Xco of 10.08 MeV. 

We then fit the 7i + ii~ mass spectrum between 3.2 and 3.6 GeV with two Gaussian- 
smeared Breit-Wigners with resolutions fixed at 10.08 and 7.83 MeV, and a second order 
polynomial background function, and with T Xc2 fixed at the PDG value of 2.00 MeV. The fit, 
shown as the curve in Fig. p], gives r Xc0 = 14.3 ±2.0 MeV, where the error is statistical. The 
fitted numbers of Xcj — ► k + h~ events are n° bs _^+^- = 720 ± 32 and n° bs = 185 ± 16, 

where the errors are statistical. 

We fit the K + K~ mass spectrum between 3.2 and 3.6 GeV to two Gaussian-smeared 
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Breit-Wigner resonance functions plus a background function that includes the possibility 
of distortions to the line shape due to /i + /i~X background events. (Because this mode is not 
used for width determination, a high precision knowledge of the mass resolution is not an 
issue.) The Xc2 width is fixed. The resulting fit, shown in Fig. ||, gives n>^_^ K+K - = 774±38 
and n°^ 2 ^ K+K - — 115 ± 13, where the errors are statistical. 

Errors in the determination of r Xc0 are caused by the uncertainty in T Xc2 , the determi- 
nation of the mass resolution, the shape of the background, the mass dependence of the 
efficiency correction, and the choice of experimental cuts. We add the estimated errors from 
these sources in quadrature and get a total relative systematic error on r Xc0 of 21%. 

Systematic errors on the branching fractions, which arise from the uncertainties in T Xc2 , 
the mass resolution, the choice of the background function, the efficiency determination, and 
the choices of the selection criteria are 11.5%, 12.6%, 12.1% and 14.7% for B(xco —* vr + 7r~), 
B(xco —> K + K~), B(xc2 — ► 7r + 7r - ) and B(xc2 — > K + K~), respectively. There are overall 
errors caused by the uncertainty of the total number of if)(2S) events and the uncertainties 
in the if)(2S) — > '"fXcj branching fractions. Adding these errors in quadrature gives total 
relative systematic errors of 14%, 15%, 15% and 16%, respectively, for B(xco ~> 
B(x c0 ^K+K-),B(xc2^K + K-)^dB(xc2^K+K-). 

For B(xcj —> Pp), sources of systematic errors include those listed above plus that 
associated with the assumption of an isotropic angular distribution for Xa,2 — » PP decays. 
Adding all of the errors in quadrature gives relative systematic errors of 33%, 67% and 56% 
for the Xco, Xci an d Xc2 states, respectively. 

In summary, we obtain the total width of the Xco to be 

T Xc0 = 14.3 ± 2.0 ± 3.0 MeV, 

where the first error is statistical and the second is systematic. The Xcj branching fraction 
results are listed in Table | and Table 0. The measured width for the Xco is consistent 
with but substantially more precise than the previous measurement 0] (the uncertainty is 
reduced from 40% to 25%). The calculations involving new factorization schemes with 
high order QCD corrections [0)0 are in good agreement with our measurement. 



Our branching fraction for Xco ~~ * PP is the first measurement for this decay, and is 
compatible with the previous upper bound |J. Our results for Xci,2 PP decays, although 
statistically limited, are consistent, within errors, with the values determined from studies 
of charmonium states formed directly in pp annihilation jnj. The calculation with mass 
correction effect Q gives much smaller value of r(x c o — > pp) than our result while the 
model considering the diquark content of the proton [||] can find result consistent with our 
measurement. 

Finally, our branching fractions for Xco,2 decays into 7r + 7r~ and K + K~ are somewhat 



lower than the existing world average [flQf| . Recent calculations of exclusive Xcj decays that 
include contributions from color-octet processes ]T5] are in generally good agreement with 
our measurements. Using our results and canceling out the common errors in the branching 
fractions, we get the ratios of the branching fractions of = 0.82 ± 0.15 and 

w Xc2 ^k \ = 1-88 ±0.51. 

B(Xc2->K+K ) 
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FIG. 3. The pp mass distribution for selected 4>{2S) — > 7pp events. 
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TABLES 



TABLE I. Branching fractions of Xcj -> _ _ . B(tp(2S) 7x c0 ) = (9.3 ± 0.8)% and 
B(tp(2S) — > 7Xc2) = (7.8 ± 0.8)% are used for branching fractions determination. 



decay mode 


jyobs 


BR(xlO~ 3 ) 


PDG @ (10- 3 ) 


XcO vr + TT - 


720 ± 32 


4.68 ±0.26 ±0.65 


7.5 ±2.1 


XcO - ^ + ^" 


774 ± 38 


5.68 ±0.35 ±0.85 


7.1 ±2.4 


Xc2 -> VT+7T" 


185 ± 16 


1.49 ±0.14 ±0.22 


1.9 ± 1.0 


Xc2 - i^+i^- 


115 ± 13 


0.79 ±0.14 ±0.13 


1.5 ± 1.1 



TABLE II. Branching ratios of \cj -> VP- B(ip{2S) -» 7x c0 ) = (9-3 ± 0.8)%, 
B(if>(2S) -► 7Xci) = (8.7±0.8)% and B(ip(2S) -> 7X c2 ) = (7.8±0.8)% are used for branching frac- 
tions determination. P* ' from this experiment, L' v ot = 0.88±0.14MeV and r*°* = 2.00±0.18MeV 

/CcU ' Acl A.c2 

are used in calculating the partial widths. 



state 


N obs 


BR(xlO" 5 ) 


T pP (keV) 


PDG [Hj BR(10~ 5 ) 


PDG @ L pp -(keV) 


XcO 


15.2 ±4.1 


15.9 ±4.3 ±5.3 


2.3 ±1.1 


< 90 




Xcl 


4.2 ±2.2 


4.2 ±2.2 ±2.8 


0.037 ± 0.032 


8.6 ± 1.2 


0.074 ± 0.009 


Xc2 


4.7 ±2.5 


5.8 ±3.1 ±3.2 


0.116 ±0.090 


10.0 ±1.0 


0.206 ± 0.022 
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